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Abstract. The generalized Floquet approach is developed to study memory effect 
on electron transport phenomena through a periodically driven single quantum dot in 
an electrode-multi-level dot-electrode nanoscale quantum device. The memory effect 
is treated using a multi-function Lorcntzian spectral density (LSD) model that mimics 
the spectral density of each electrode in terms of multiple Lorentzian functions. For 
the symmetric single- function LSD model involving a single-level dot, the underlying 
single-particle propagator is shown to be related to a 2 x 2 effective time-dependent 
Hamiltonian that includes both the periodic external field and the electrode memory 
effect. By invoking the generalized Van Vleck (GVV) nearly degenerate perturbation 
theory, an analytical Tien- Gordon-like expression is derived for arbitrary order multi- 
photon resonance d.c. tunneling current. Numerically converged simulations and 
the GVV analytical results are in good agreement, revealing the origin of multi- 
photon coherent destruction of tunneling and accounting for the suppression of the 
staircase jumps of d.c. current due to the memory effect. Specially, a novel blockade 
phenomenon is observed, showing distinctive oscillations in the field-induced current 
in the large bias voltage limit. 
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1. Introduction 

Electron transport of a quantum system in the presence of time-dependent external 
fields has been studied by various approaches, often leading to observation of new 
phenomena and applications. An early experiment conducted by Dayem and Martin in 
1962 p] studied photon-assisted tunneling (PAT) processes in superconductor-insulator- 
superconductor hybrid structures, in which temporally periodic fields were applied to 
the source and drain. Subsequently, Tien and Gordon proposed a theoretical model of 
the PAT in 1963, suggesting that a time-dependent periodic external field can produce 
distinctive sideband structures of the electron density in the source and drain [2J. 
Particularly, in the past two decades, due to the advent of nanotechnology, the effect 
of time-dependent fields on the electron tunneling through nanoscale devices has been 
extensively investigated both experimentally and theoretically. On the experimental 
side, PAT has been studied in various nanoscale systems, including GaAs/AlGaAs 
quantum dots [31 E] and single-donor quantum dots in semiconductor nanostructure 
[5]. Especially, the staircase d.c. current as a function of bias voltage (the I — V 
characteristics) in the Coulomb blockade regime [H] and the current oscillation around 
zero bias voltage due to external field [7] have been observed. On the theoretical side, 
various treatments of PAT for nanoscale devices have been formulated based especially 
on quantum master equation approaches [HJ E] and the scattering theory in the context 
of the non-equilibrium Green's function (NEGF) method [101 [TlJ [121 EES E]> including 
multiple photon assisted tunneling phenomena [T5\ US] and Non-Markovian memory 
effect [TTJ . The PAT absorption/emission sideband structures revealed in these studies 
agreed with the experimental findings. 

The Floquet theory has undergone extensive development and generalization in 
the last few decades [18], including the Floquet matrix method [191 1201 [21], many- 
mode Floquet theorem [22], and Floquet-Liouville super-matrix formalism [23] . More 
recently, it has been extended to study quantum interference of periodically driven 
superconducting qubits [23]. In addition, the generalized Van Vleck (GVV) nearly 
degenerate perturbation theory can be used to analytically study near-resonant multi- 
photon processes for few level systems J2J3 [2S] . In the past decade, the Floquet approach, 
within the wide-band limit, has also been widely adopted in the study of periodically 
driven electron transport processes involving nanoscale quantum devices [27]. 

In the present work, we extend a recently formulated generalized Floquet theory, 
amenable to the electrode memory effect, to treat time-dependent electron transport 
phenomena through a periodically driven single quantum dot. A nanoscale quantum 
system may contain an externally driven central junction, which usually are quantum 
dots [271 [28] or single molecules [291 [301 [31] . Theoretically, nanoscale quantum devices 
are commonly treated in terms of tight-binding models. In the adiabatic limit, each 
quantum dot in the tight-binding model is usually endowed with a single electronic 
level that allows one electron at a time to tunnel through. In the general non-adiabatic 
situation, each quantum dot may possess multiple coupled electronic levels (as well 
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Figure 1. Energy diagram of multi- level tight-binding model of a single quantum dot 
in contact with two electrodes (left and right) . The dot possesses discrete energy levels 
£i in the presence of an external field A(t). The electrode are free electron Fermi gases 
with the electrochemical potential /i£ and /ir respectively. The coupling between the 
SQD and the electrodes are represented by the spectral density f ^(e) and f n(e). 



as vibrational ones for molecules junction), thus permitting multiple electron tunneling 
pathways that can lead to a variety of constructive and destructive interference patterns. 
In addition, the presence of time-dependent external fields applied directly to the central 
junction can open up new pathways due to single- and multiple-photon resonance 
processes. The generalized Floquet approach in this paper is formulated for a single 
quantum dot junction endowed with N noninteracting electronic states [12] as shown 
in Fig. [TJ This model does not permit direct electron transition between dot levels, 
however each level still possesses a non-zero dipole moment. The spectral density of the 
electrode-junction coupling T L / R (e) is represented by the Lorentzian spectral density 
(LSD) model in terms of a sum of multiple Lorentzian functions. The same LSD model 
has also been adopted in the generalized Floquet approach developed for the electrode- 
multi-dots-electrode quantum system [32J . 

The remaining parts of the paper are arranged as follows. In Sec. 2, we describe the 
general formulation and the generalized Floquet approach for the electrode-multi-level- 
electrode system driven by a time-periodic field. In Sec. 3, we consider the special case 
of the symmetric single-function LSD (SS-LSD) model and a single-level quantum dot. 
The memory effect on the multi-photon resonance PAT is studied using the GVV nearly 
degenerate perturbation theory. In Sec. 4, we compare the analytical GVV results and 
the corresponding exact Floquet calculation, especially in the limit of large bias voltage. 
Finally, a summary is given in Sec. 5. 

2. General Formulation 

2.1. Driven multi-level model and scattering formalism 

Consider a single quantum dot (SQD) in direct contact with a left electrode (source) 
and a right electrode (drain), as shown in Fig.l. The Hamiltonian of the electrode-SQD- 
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electrode system can be written as [HlJ EZ], 

H(t) = H c (t)+ H e + H', (1) 

l=L,R 

where the external field is applied only on the SQD. Specifically, we consider a 
periodically driven single quantum dot that possesses N non-interacting electronic 
states. In the tight-binding model, the Hamiltonian of the central quantum dot takes 
on the expression 

N N 

H c {t) = e i d Ui + A (t) E Vi4 <k, (2) 

i=l i=l 

where di(d\) and \i\ are, respectively, the energy, the annihilation (creation) operator 
and the dipole moment of the i-th electronic state |i) [6j[33]. The external field A(t) is 
a periodic function of time, i.e., A(t + T) = A(t). The left and right electrodes are free 
electron Fermi gases (in thermal equilibrium) and can be described by the Hamiltonians 

H e = Y e tq4 q c eqi ? = L,R, (3) 

where c^ q (c\ ) is the annihilation (creation) operator for the electron states \£q) 
associated with the energy ee q . The contact Hamiltonian between the SQD and two 
electrodes is given as 

H ' = E ( E V ^ c lA + E V mAA + h.c. (4) 

i=l \ Lq Rq / 

where Vi q ,i and Vr^^ are, respectively, the corresponding coupling parameters. 

Within the framework of the non-equilibrium Green's functions (NEGF) method, 
the corresponding single-particle retarded Green's function can be computed via the 
relation 

1 r°° 

G(t, e) = G(t + T, e) = — / e tet ' /R U(t, t - t')dt', t > t' > 0, (5) 

Jo 

where U(t + T, t + T) = U (t, t ) is a N x N time-dependent matrix, representing the 
underlying single-particle propagator for the SQD. The fcth-order Fourier coefficient of 
the retarded Green's function G(t, e) is given as 

G^(e) = ^£G(t,e)e^ t dt, (6) 

which denotes the single-particle scattering associated with the absorption or emission 
of | A; | photons at the incident energy e [27]. The transmission coefficient of the electron 
tunneling for the |/c|-photon process is [TU] 

1 N 

1%rl^) = 4 E f + *MlWe) |(i|G ( %)|j>f , (7) 

y=i 
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where r^i(e) is the spectral density for the i-th electronic state. The total transmission 
coefficient of electron tunneling through the SQD is 

oo 

T LR/RL {e) = T lIrl(£) (8) 

k=— oo 

The time-ensemble averaged current (i.e. the d.c. current) can be expressed in terms of 
the total transmission coefficients as follows 

- e 



h 



{T LR (e)f R (e)-T RL (e)f L (e)}de, (9) 



oo 



where = 1/{1 + e < - e ~ ML < /i? ^ fcs7 ~} are the corresponding Fermi-Dirac function at 

the temperature T describing the electron energy distributions in the electrodes L and 
R, with Hl/r being the respective electrochemical potentials. The bias voltage across 
the SQD is the difference V = fi R — hl- In the zero temperature limit, i.e. T = 0, the 
Fermi-Dirac function reduces to a step function bounded by the bias voltage window 
[Vl,V R ]- 

2.2. Lorentzian spectral density model 

For a non-interacting iV-level SQD, assuming the electrodes are initially in thermal 
equilibrium [32], the corresponding single-particle propagator U(t,t Q ) is composed of N 
independent components [/;(£, t ) = (i\U(t, t )\i), i = 1,...,N, governed by the integro- 
differential equation 

ifrjUi{t, t ) = [a + mA(t)]Ui(t, h)-\Y. I - W, t )dt', (io) 



L,R 



where Ti^t — t') is the response function (memory kernel) that can be expressed as 

2 

h 



q 

in terms of the coupling parameters Vt q ^. Here, we have taken into account the two spin 
states of electrons. The spectral density l\i(e) = ^Y^i{t)e iet ' n dt can be written as 

f, ii (e) = ^47r|^ ii | 2 (5(e-e, (? ), (12) 
which a collection of delta functions at individual energy levels weighted by the 
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corresponding coupling parameters \Vi q ^ 

In the Lorentzian spectral density (LSD) model, the spectral density l\i(e) is 
considered as a linear combination of M Lorentzian functions, i.e., 

M a e U M (b e ) 2 

( £ - <4j + (%) t^l ( £ ~ ^k) + ( & fk) 

where af k , bf k and of k are the fitting parameters for the k-th Lorentzian function. The 
spectral density l\i(e) is usually a smooth function of energy e and can be readily 
mimicked by a finite number of Lorentzian functions. Here Tf k (= af k /6f k ), 6f k and of k 
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denote the coupling strength, the band-width and peak position, respectively, of the 
k-th Lorentzian function for the i-th state. From Eqs. ( 11 )-( 13 ), the response function 
can be explicitly expressed as 



M 

— Y 



2nh 



M 



k=l 



a 1 b l 
"ik°ik 



l \2 



e(t-t')/H de 



Z> 2h aikC 



(14) 



k=l 



In the wide-band limit (WBL), bf k —> oo, thus F^^e) —> ^ k rf k = is independent of 
energy e. As a result, we have F^(t — if) — x S(t — t'), which is free of the memory 
effect. 

By introducing the auxiliary functions 

f 



a[ k exp[-2( i a ik 



ibi)(t-t')/H] x Utf ,t Q )dt' , (15) 



for i = L, R, i = 1, N and k = 1, M, Eq. (10) can be recast as a set of 2M + 1 
coupled ordinary differential equations of L/j(t, t ), YJ L (t,to) and Y^it,^), i-e., 



ih 



d 
dt 



( Ui(t,t ) ^ 




( Ui(t,t ) \ 


Y^(t,t ) 


= Hi{t) 








\Yi R (t,t ) J 



(16) 



subject to the initial conditions Ui(t ,t ) = 1, YJ L (t ,to) = 0, and Y| R (to,to) = 0. Here, 
the effective Hamiltonian "Hi(t) take on the form 

( e x + ^A{t) Qt Q R \ 

Qf f Ef - iB? 



where 



and 



\ 



Q 



I 



or 



o \ 



\ 



(17) 



(18) 



(19) 



J iM 



J 



The effective Hamiltonian is a time-periodic non-Hermitian matrix. Eq. (16) may 
effectively be considered as the governing equation of a periodically driven N 
noninteracting levels coupled to M unstable levels modeling the spectral densities of 
states in the left and right electrodes. The individual level index i does not play a role 
in the following formulation and, therefore, will be dropped for simplicity. 
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2.3. Generalized Floquet approach 



By invoking the generalized Floquet theory, the solution of Eq. (16) can be written 

as 

y(t) = $(£) x e ~ iAt/h , (20) 

where A is a diagonal matrix composed of complex quasienergies X a and is a time- 
dependent periodic matrix, $(t + T) = composed of the corresponding Floquet 
quasi-states 4> a (t), 4> a (t + T) = (fi a (t), satisfying the eigenvalue equation, 

H{t)-iH^)(f> a (t) = \ a <f> a (t). (21) 



dt / 

in an extended Hilbert space [HI [20]. By expanding the Floquet states <f> a (t) 



2M+1 oo 



= E E A>> xe ^ ( 22 ) 



=1 n=— oo 



Eq. (21) can be recast as a time-independent quasienergy equation 

2A/+1 oo 

E {(o^W-^ + mMMj^^O. (23) 

/3=1 ?i=— oo 

and its concomitant adjoint equation 

2A/+1 oo 

£ J2 {H?* t(m ~ n) |/?} - (A^ + mhw)S afl 5 mn } $g = 0, (24) 

/3=1 n=— oo 

for a = 1, 2M + 1 where |a) and \(3) indicate the orthogonal unperturbed eigenstates 
with A(t) = and U {n) = ± H{t)e inult dt. 

The fundamental solution W(t,t ) = ^(t)^ -1 ^) = $(t)e -iA ( t - t|, )$- 1 (to) can be 
obtained by solving Eq. (21). Given the initial condition, Ui(to,to) = 1, the desired 
single-particle propagator can be obtained Ui(t,t ) = Uu(t,t ). In practical, we can 
solve the quasienergy equation numerically by truncating the Floquet Hamiltonian. 



3. Symmetric single-function LSD (SS-LSD) model and a single-level 
quantum dot 

In this section, we consider the symmetric single-function LSD model and a single- 
level quantum dot. In particular, GVV nearly-degenerate perturbation theory is adopted 
to derive a Tien-Gordon-like expression for the d.c. current at single- and multi- 
photon resonance conditions. The resultant derivations can be readily extended to 
non-interacting multiple level cases. 
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(a) (b) 




Figure 2. Energy diagram of the SS-LSD model and a driven single-level quantum dot 
are shown in (a) , where the shade areas represent the single Lorentzian spectral density. 
The effective three-level system (b) are consist of effective levels (thick dash lines), 
representing the center of Lorentzian function: a L = a R = a. The rt-photon resonance 
condition and m-photon absorption inelastic scattering process are also shown. 



3. 1 . General formulation 

We assume that the single dot level has the energy eo and dipole moment /j,q = 1. 
The SS-LSD model depicts the left and right spectral densities with the same Lorentzian 
function, i.e., 

r i(e) = r R (e) = f ( C ) = ^ = (25) 



where T = a/b. Within the SS-LSD model, Eq. (16) is reduced to 

d_( U(t,to)\ (eo + A(t) iyfi/2\( U(t,t ) 
m dt [ Y+(t,t ) ) [ -i^/a~j2 a -ib )[ Y+(t,t ) 1 

where Y ± (t,t ) = (Y R (t,t ) ± Y L (t, t Q ))/y/2 and A(t) = Acosut. Noting that 
Y~(t, t ) = Vt > t since it is decoupled from both U(t, t ) and Y + (t, t ) and possesses 
the initial condition Y~(to,to) = 0. Consequently, the transmission coefficient and the 
d.c. current can be written, respectively, as 

T RL (e)=T LR (e) = T(e) = i^f (e + ktiw)t(e)\G^(e)\ 2 (27) 

k 



T(e)[f R (e) - f L (e))de. (28) 



and 

'-5 

It is seen that in the WBL the transmission coefficient T(e) can be further reduced to 
T WBL {e) = \v^\G {k \e)\\ (29) 



for r(e) ->• T = a/b as b -»■ oo, cf. Eq. (25). 
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Near the n-photon resonance condition a — eo ~ nhu and introducing the rotating 
frame 



U{t,t ) 
Y+(t,t ) 



K(t) 



X'(t,t ) 
Y'(t,t ) 



where 



n{t) 



-in(t) 





-incut 



(30) 



(31) 



with fi(t) = | J A(t')dt' = -^-sinut, Eq. (26) can be transformed to 



I X>(t,t Q ) 
dt I Y'(t,t ) 



X'(t,t ) 
Y'{t,t ) 



where 



and 



V(t) 



e 

a — nhu — ih 



_j e -i(Q(t)~nujt) Q 



(32) 



(33) 



(34) 



y( m ) 



Here, £ = \/a/2 is a small parameter associated with the weak coupling and finite band 
width. Consequently, the single-particle propagator can be computed using the relation 

U(t,t') = [Jl(t)U'(t,t')rt(t')] n . (35) 

Invoking the relation e ixsmult = Yl m Jm{x)e' imuJt , x = A/hu and J m (x) is the m-th order 
Bessel function of the first kind, V(t) can be Fourier expanded as V(t) = J2 m V^ m 'e' 

it/ n _ m (x) 

At the n-photon resonance condition, the off-diagonal part of the rotated Hamiltonian 
contain the couplings J n ± m {x) corresponding to the |m|-photon absorption (emission), 
as shown in Fig. [2j 

3.2. Generalized Van Vleck (GVV) nearly- degenerate perturbation theory 

Near the n-photon resonance condition a — e ~ nhu, with the aid of the generalized 



-irriLdt 



(36) 



Van Vleck (GVV) perturbation method (see Appendix), Eq. (32) can be further 
approximated as 



ih 



d / x%t Q ) 



n 



Of \ Y%t ) J '" GVV \ Y'(t,t ) 
where the nth-order resonance GVV Hamiltonian takes on the form 

eo - i 2 $n i£J n {x) 



X'(t,t ) 



H 



GVV 



-i£J n {x) e -ib + £, 2 5 r , 



(37) 



(38) 
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in which the off-diagonal terms ±i£,J n (x) are the first-order correction with respect to 
the small parameter £ while the diagonal terms are the second order correction with 

i E q ^4*W J ™& if n = 

m^O m^O 

It can be seen that the real part of the second order correction £ 2 5 n is responsible for the 
level shift [24] and the imaginary part, together with the first order correction i£J n (x), 
is responsible for the narrowing of the PAT sidebands. 

As shown in Appendix, Eq. (A18), the underlying retarded Green's function can 
be written as 

Jm\ x )Jm+k\ x ) (40) 




e — e — mhu + A n + iE Ti 



if n = 

I E (dg^nVJ*) if n > 1 ( 41 ) 



and the corresponding width is 



1 * ^ — \ b 2 2 / \ r a , , 

2 ^ (mM 2 + & 2 +m 2 = 26' 

m=— oo v 



Furthermore, by substituting Eq. (|40j) into Eq. (|8j), we can derive an analytical 
expression for the transmission coefficient 

m=— oo v 

where the tunneling rate of the m-th PAT sideband is 

00 °° b 2 T 2 (r) 

i m {e) = £ r(e + = r ]T (e + fc ^;^ 2 - (44) 

It is seen that in the WBL 7 m (e) — >■ T, A n — >• and E n — > |), leading to 

j 00 p2 

7W(e) = - ^ J 2 t (x)x- p2 . (45) 

4 m t"oo (e - eo - ™M 2 + T 

Both TvysL(e) and Tcw(e) are composed of a collection of sidebands. However, the 
GVV sidebands exhibit additional shifts (— A n ) and are also narrower in width (S„ < |) 
due to the memory effect. 

From Eqs. (28) and (43), a Tien-Gordon-like formula for the corresponding d.c. 
current can be written as 

law = Yl J ™^) x 7 ™' ( 46 ) 



m=— 00 
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Figure 3. The amplitude dependence for (a) the level shifting A„ and (b) the sideband 
width S„ are plotted with respect to various resonance conditions n = (a — eo)/huj 
for the SS-LSD model with a = and b = 1 Huj . The thick gray line indicate the 
corresponding values in the wide-band limit (i.e. b — > oo). 



which is composed of an infinite number of weighted contributing partial currents 

J m = ^ r 7 T{ T m }i , 2 ^ 2 x [f R (e) - f L (e)]de, (47) 

for m = — oo, • • • , +oo. For the iV-level quantum dot, the total d.c. current can be 
expressed as a summation of individual currents through different states. 



4. Results and discussions 



In this section, we present and discuss the memory effect on the electron transport 
processes of the electrode-SQD-electrode device in the presence of a periodical field, 
based on both the numerically converged results and approximate GVV results for the 
special case of symmetric single-function LSD. First, we study the level shifting and 
width narrowing of the PAT sidebands as a function of the external field amplitude. 
Second, we compare the transmission coefficient for the SS-LSD plus a single-level 
quantum dot and that in the wide-band limit (WBL). Third, we study the memory 
effect on the staircase jumping of d.c. current. Finally, we show the field-induced 
current oscillations as a function of the gate voltage (i.e. e ) at the large bias voltage. 
In our calculations, the frequency of the driving field {A(t) = A cos out) is fixed at 
hw = 10 eV, which is also used as the unit of energy, and the field amplitude A is 
allowed to take on arbitrary values. The electrode-single dot coupling strength is chosen 
as T L = T R = T = 0.1 hw. 



Figure |3| shows the results based on Eqs. (41) and (42), depicting the memory 
effect on the level shifting A n and width narrowing S n , as a function of the amplitude 
x (x = A/hu), at the n = 0, 1, 2, 3-photon resonance conditions. It is found that in 
the small amplitude limit i < 1, both A n and S n are dominated by the term |Jo(x)| 2 , 
corresponding to m = —n in Eqs. (41) and (42). At x — 0, the level shifting and width 
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Figure 4. The contour plots of the transmission coefficients T(e) are plotted as a 
function of the incident energy e and the external field amplitude x = A/Huj for the 
WBL ((a) and (b)) and the SS-LSD model with b = 2 Hlj ((c) and (d)). We compare the 
WBL and SS-LSD results for the resonance conditions n — and n — 2. The inserts 
show the corresponding energy diagrams. The vertical dash lines indicate the roots of 
Jq(x) (left) and J2(x) (right) that are related to coherent destruction of tunneling. 



narrowing can be computed as 



T —nbhhj 



and 

_ r b 2 



> for all n 



n 2 (nhw) 2 + b 2 

respectively. In the large amplitude limit x ^> 1, \ J m (x)\ 2 ~ (2/nx) cos 2 [x — (2m+l)7r/4], 
thus A n (n 7^ 0) and S n are inversely proportional to x - they both become smaller as 
x becomes larger, with the former goes to zero much faster because of the cancellations 
between opposite signs arising from the summation index m in Eq. (41). The observed 
behaviors of the level shifting and width narrowing, as shown in Fig. [3] are important 
to understand the properties of the electron transport transmission coefficients and the 
resultant d.c. currents that are to be discussed below. 

Depicted in Fig. [4] are contour plots of the electron transport transmission 
coefficients, as a function of the incident electron energy e and the field amplitude 
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e/hw 



e/tuo 



Figure 5. The numerical (green dashed curves) and GVV (red solid curves) results 
of transmission coefficients are shown as a function of the incident energy e for the 
SS-LSD model (6 = 1 Hw) with respect to various resonance conditions (a) n = 0, 

(b) n = 1, (c) n = 2. The corresponding energy diagrams are as the inserts and the 
external field amplitude is x = 2.40. We label each peak by the m index , including the 
vanishing peaks at m = due to CDT. The scaled insets on the right corner in (b) and 

(c) indicate the level shifting near e = 0. The corresponding transmission coefficients 
for the WBL are shown in the right panels. 



x for the SS-LSD parameters o = (spectral peak center) and b = 2 hw (spectral 
width). The calculations were done based on Eqs. (23) and (24). It is found that the 
PAT sideband structure contains patterns of narrow peaks at the multiple integers of 
the applied field frequency hw, i.e., e = e + khw, k = 0,±1,±2, corresponding 
to various multi-photon absorption and emission processes. Furthermore, in the WBL 
(Figs. [4j^a) and (b)), the transmission coefficient contours display the same sideband 
structures that are centered at different resonance energies, e = hw (n = 0) and 
e = 2hw (n = 2), respectively, cf. Eq. (45). However, in the SS-LSD model (Figs. 
|l]^c) and (d)), the transmission coefficients reveal a very different sideband patterns for 
different incident energies - here the SS-LSD sidebands are symmetric for n = but 
asymmetric for n — 2, showing the results of the level shifting and width narrowing 
because of the memory effect, cf. Eq. (43). 



Figure |5| further shows numerical and GVV (green-dashed and red-solid, 
respectively, in the left panels) as well as the WBL calculations (the right panels) of the 
non-overlapping peaks (m = 0, ±1, ±2, ■ ■ •) in the transmission coefficients as a function 
of the incident energy e for n = 0, 1, 2. It is found that the PAT sideband is completely 
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Figure 6. The d.c. current / is plotted as a function of bias voltage V — [Ir — fi^ 
with fixed (II — in the upper panel and its corresponding transmission coefficient is 
shown in the lower panel. The results for the SS-LSD model (eo = cr = and 6 = 1 fiu?) 
is compared with the WL results (gray dash curves) at x = 2 Twj. 



suppressed at m = (i.e., e = —nhu), a manifestation of the CDT phenomenon 152]. 
The inserts on the upper right corners in the panels (b) and (c) display the enlarged 
peaks near e = 0, revealing significant memory effect. Clearly, as shown in the left 
panels of Fig. [5j each n-photon [n = 0, 1,2) resonance PAT transmission coefficient is 
composed of a sequence of non-overlapping sidebands (m = 0, ±1, ±2, ■ ■ ■), respectively, 
corresponding to e/hw = —n + m with the amplitudes approximated as 

, , r 2 Jl(x) 

7(n,m) x x 



where 



4(A2 +S 2) (- n + m )2 + i> 

7(n, m) = > t — -. 

v ' ^ i-n + m + k) 2 + 1 

k=—oo 

The n — 1, 2-photon PAT sidebands in Figs. (5^b) and are not symmetric about 
m = due to the asymmetric weights of j(n, m) and 7(n, — m). This is in contrast to 
the always symmetric sideband structures about m = in the WBL, as shown in the 
right panel in Fig. [5j 

In Fig. |6j the upper panel shows the staircase feature of the d.c. current J as a 
function of the bias voltage V = [Mr — Ht (here fiL = 0), whereas the lower panel shows 
the sideband structure of the underlying transmission coefficient T(e) as a function of 
the incident energy e. The d.c. current staircase is a result of the increasing number of 
the sidebands that are located within the bias voltage window [//£,//#]. In general, as 
the bias voltage window becomes bigger, there are more non-vanishing PAT sidebands 
contributing to the d.c. current in the WBL than in the SS-LSD model. In the 
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Figure 7. The d.c. currents 7 and the field-induced current T — To are plotted as a 
function of the level energy (eq) for a large bias voltage = — /!/? = 10 [fiw]) in the 
SS-LSD model (a = and 6=1 shown in (a) and (b). The insets shows the 
corresponding energy diagrams. We also plot the d.c. current in the WBL for the zero 
bias voltage case (fi^ = ~ A'-R = 0.1 Hlu) in the panel (c). 



latter, the memory effect quickly suppresses the higher |m| PAT sidebands, leading 
to a much smaller d.c. current. The insert in the upper corner illustrates how the SS- 
LSD (6=1 fajj) staircase feature (solid curve) gets smoothed out in the WBL (dashed 
curve) . 

Figure [7] shows the memory effect in the SS-LSD model on the electron tunneling 
blockade phenomenon in the large bias voltage limit. Specifically, the d.c. current I 
calculated using a bias voltage window [p, L = 10 hu,fiR = —10 Hlu] is shown as a 
function of the reference gate voltage eo for different amplitudes A of the external field 
(with a = and b = 1 hu fixed in the calculations). Fig. [7^a) shows that (1) the d.c. 
current feature (as a function of e ) strongly depends on the value of x{= AjfoJ) and (2) 
the d.c. current drops to zero as soon as the reference gate voltage eo moves outside the 
SS-LSD width b. The latter finding is responsible for the d.c. current blockade caused 
by the memory effect in the SS-LSD model. The amplitude dependent behavior seen 
in Fig. [Tj(b) is complete absent in the WBL where the d.c. current at very large bias 
voltage can be computed explicitly as 
- e irT 

WBL = XT' 

which depends only on the coupling strength T (here Iwbl = 5.772 x 10~ 2 [e/h] for 
T = OAhu). By taking the memory effect into account, the tunneling current is 
effectively blocked between the SQD and the electrode because of the limited accessible 
energy bandwidth in the SS-LSD model. The tunneling coupling between the SQD 
and the electrodes is strongest when the energy level eo coincides with a, see the inset 
in Fig. [7|a). In addition, Fig. ^b) demonstrates that there exist significant oscillatory 



16 



Coulomb blockade features in the field-induced current I — Iq where Jo is the d.c. current 
in the absence of external fields. The large bias voltage limit considered in Fig. [7Fa) 
and [7](b) clearly reveals a very different mechanism (due to the memory effect) for the 
d.c. current Coulomb blockade phenomena from the mechanism responsible for the 
small bias situation in the WBL, in which the electron can tunnel through the SQD 
only when one of the sideband resides inside the small bias voltage window, as depicted 
in Fig. [7Fc). The interplay of the band width b and the bias voltage window [hl-iI^r] 
may result in the enhancement or suppression of the d.c. current by manipulating the 
gate voltage e . 

5. Conclusions 

In summary, we have developed a generalized Floquet approach, including the 
memory effect, for the treatment of electron transport process on a periodically-driven 
single quantum dot system with multiple noninteracting levels. Of particular interest, 
we have considered the symmetric single-function Lorentzian spectral density (SS- 
LSD) model for the electrodes and derived analytical expressions for the transmission 
coefficient and the d.c. current under the multi-photon resonance condition by the 
generalized Van Vleck (GVV) nearly-degenerate perturbation theory. The Tien-Gordon 
formula has been extended to include the memory effect and the multi-photon resonance 
processes, in particular, resulting in an effective multi-level model of a single quantum 
dot (one single level with multiple sidebands). The memory effect on the transmission 
coefficient and the d.c. current has been analyzed at the nearly-degenerate resonance 
conditions and for different external field amplitudes. Numerical simulations of the 
transmission coefficients have shown that some multi-photon PAT sidebands can be 
suppressed by the memory effect. We have also shown that the memory effect on the 
staircase feature of the d.c. current is closely related to the sideband width narrowing. 
Furthermore, It has been observed that the electron tunneling may be blocked for large 
bias voltage case if the gate voltage is moved outside the band width of the electrode 
spectral density function. The field-induced current oscillation may be produced and 
manipulated by applying periodic external fields, thus, enabling the enhancement or 
suppression of the d.c. current at certain gate voltages. 
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Appendix. Derivations of Eq. (38) and the corresponding single-particle 
propagator 



The time-dependent equation, Eq. (32), can be written in terms of the Floquet 
quasi-states \an) = \a) x e mujt in extended Hilbert space, leading to the Floquet 
Hamiltonian 



where 
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where I being a 2 x 2 identity submatrix and H and V are, respectively, composed of 
2x2 sub-blocks shown in Eq. (33) and Eq. (36). 

By invoking the nearly degenerate GVV perturbation theory, an effective n-photon 
resonant 2x2 Hamiltonian Hew and its eigenstates can be expressed as 



n 



GVV 



52 Cn 



(m) 



m=0 

oo 



i^> = E^>i m) > 



m=0 



In terms of extended Hilbert space, the zeroth order eigenstates are 



|^ 0) ) = |1,0), |^> = |2,0>. 



(o)\ 



The first-order perturbation terms are, respectively, 





-iJJx) 



iJ n (x) 





(A4) 
(A5) 

(A6) 
(A7) 
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and, 

I (1)\ ~^Jn+m\X) , . | (l). iJ n+m (x) . . / a q \ 

The second order perturbation terms for %gvv can De written as 

S n 

5 n 

where 



^ (2) = ( "5 r ° I , (A9) 



■^■f — z6 + mtko 

< E WW^Jl^) if « = ( A1 °) 

E (m7fcj) 2 +6 2 ^n+m( X ) + * E (^Jp+p ^n+m ( X ) ^ n 7^ 

Hence, the GVV effective Hamiltonian, up to the second-order, can then take on the 
form 



now* i t0 ,; 7 \:- . ^rrL I, (ah) 



£0 - £ 2 $n i£Jn{x) 

-i£J n (x) e -ib + £ 2 5 n 

compared to the corresponding effective Hamiltonian in the rotating wave approximation 
(RWA), 



W»n I a | ./° , ' tiUX l ]. (A12) 
-i£J n (x) t Q -ib j 

In the weak coupling limit and with a finite band width, the eigenvalues of the 
GVV effective Hamiltonian T-Lgvv can be approximated as 



~ e — A n — zc n 
A_ ~ cr — nfta> — z6 + A„ + iS r 



(A13) 



where 



if n = 

A " " i IE (^Tp4 2 +m (x) if n > 1 ( A14 ) 



and 



r °° o 2 

S "=2 ( m ^) 2 + 6 2J " +m(a;) ' (A15) 

m=— oo v ' 

Hence, we can derive the fundamental solution by U(t,t') « lZ(t)UGvv(t, t')7V(t') and 
the single-particle propagator U(t,t'), cf. Eqs. (30) and (35), 

(7(M') « e -*n(*) e -*Mt-f)/R e in(*0 (A16) 

In conjunction of e lxsmw * = J m (x)e* mw *, we derive the Green's function, cf. Eq. (pi), 



G(M)K £ £ ^M xe -,. (A17) 

-oo m=— oo 
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and its k-th order Fourier coefficients, cf. Eq. rt6 



are 




G«(c)« £ 



e — e — m/iu; + A n + iH 



(A18) 



m=— oo 
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